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DETERMINATION OF THE THEORETICAL PRESSURE DISTRIBUTION FOR 

TWENTY AIRFOILS 
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SUMMARY 

This report gives the theoretical distribution of pressure 
at lift coefficients of 0, 0.5, 1.0, and 1.5 for SO airfoils, 
calculated on the basis of a rigorous potential theory of 
arbitrary airfoils. It also provides tables from which 
the characteristics of the airfoils for any angle of attack 
in dimensional potential flow are readily calculable. 
The theoretical values of the angles of zero lift, the lift and 
moment coefficients, and the ideal angles of attack are 
listed and some comparisons with experiment are indi- 
cated. Some of the well-known characteristics and 
properties of airfoils are accounted for in terms of the 
theoretical pressure-distribution curves. Qualiiatixe de- 
ductions are made concerning the causes of breakdown 
of potential flow and the efficiency of the airfoil in viscous 
flow. The results presented may be of value in pre- 
dicting structural loads and also in a correlation of theo- 
retical pressure gradients with profile resistance. 

INTRODUCTION 

Until recently the theoretical distribution of pressure 
around airfoils could be determined only for the so- 
called "theoretical" airfoils. Indeed, only in the par- 
ticular case of the Joukowsky airfoils is the calculation 
not unduly laborious. (See references 1 and 2.) 
The theoretical airfoils, which are defined by special 
mathematical transformations, have, however, seldom 
been employed in practice. Their use in a precise 
study of pressure distribution has in fact been due more 
to necessity than to desire. The distribution of pres- 
sure for mathematically "thin" airfoils (i.e., the airfoil 
is represented by the mean-camber line) can be ob- 
tained, at least approximately, by the processes given 
by Munk, Glauert, and Theodorsen (references, 3, 4, 
and 5). In another report (reference 6) Theodorsen 
developed a theory readily applicable to arbitrary 
airfoils. This theory was extended by Theodorsen 
and Garrick in reference 7, in a report which gives a 
unified treatment of the 2-dimensional potential flow 
around airfoils of any shape. The treatments given 
in references 6 and 7 avoid approximations in the anal- 



ysis, and are referred to for details of the underlying 
theory of the results of the present paper. 

The differences exhibited by airfoils in potential 
flow, as well as the differences between the actual and 
ideal cases for a particular airfoil, can, of course, be 
critically studied only if the ideal case is known. 
Furthermore, it is only on this basis that the assump- 
tions of the theory itself can be critically analyzed 
and modified. It is therefore believed that an existing 
gap in aerodynamical literature will be, to some 
extent, bridged by the publishing in the present paper 
of convenient tables and curves of the theoretical 
results for a number of co mm only used and related 
airfoils. 

A knowledge of the theoretical distribution of pres- 
sure for an airfoil is, undoubtedly, a major factor in 
making it ultimately possible to predict accurately the 
behavior and efficiency of the airfoil under actual con- 
ditions, for the theoretical changes along the surface 
from pressure to velocity and from velocity to pressure 
are very significant in the determination of the drag 
characteristics. Knowledge of the theoretical results 
is of considerable value, too, for guiding experi- 
mental work whenever the measurements are rather 
critical, and such information also directs attention to 
the significance and interpretation of differences 
between theory and experiment. 

Unfortunately, because of lack of sufficient accurate 
experimental data, comparison cannot be made directly 
with wind-tunnel results except in a few cases. In 
reference 7 an interesting comparison was given between 
theory and experiment of the pressure distribution 
around the N.A.C.A.-M6 airfoil at 12 different angles 
of attack. Reference 8 may be referred to for quali- 
tative experimental results for five additional airfoils. 
A more accurate experimental study of pressure dis- 
tributions is in progress at the present time at the 
N.A.C.A. laboratories. 

A part of the following work was undertaken at the 
request of the Bureau of Aeronautics, Navy Depart- 
ment, for use in work on structural loads. 
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In making the calculations the author was ably 
assisted by Miss Alyce V. Rudeen, of the Committee's 
staff. 

SDMMABT OF FORMULAS USED 

The formulas used to obtain the results presented in 
the tables and curves are developed in references 6 and 
7. A sample calculation for the N.A.CA.-M6 airfoil 
with a comparison with experimental results, as well 
as explanatory figures and diagrams illustrating the 
use of the formulas, is given in reference 7. The 
following list presents the symbols employed and their 
definitions: 

SYMBOL DEFINITION 



(1) 

(2) 8 

(3) # 



y) See discussion of the choice of axes in a 
following paragraph. 



(4) * 

(5) *' 

(6) V 

(7) <t> 

(8) to 

(9) « 

(10) (3 

(11) k 

(12) 4 



2 sin 2 5=^+-v^+^ 
where p = l — (| 



2 sinh 2 \fr= —p + -yjp i +y 1 . Since y is gen- 
erally small for airfoils, the following 
equation may be preferable: 

sinh^= 2 Near the leading (or 

trailing) edge ip is given approxi- 
mately by 4> = where a- is the radius 

of curvature at the leading (or trail- 
ing l ) edge. 

^ — h{ mcoti ^~ d4> 

See appendix 3 of reference 7 for method 

of evaluation. 
Obtained graphically from the e, 0 curve. 

(Denoted in reference 7.) 

Obtained graphically from' the f, 0 curve. 

(Denoted ^ in reference 7.) 

4>=e+e. 

A constant: to — y uV (<£) dci 

^7T n 



.Angle of attack with respect to the x axis. 
The angle of zero lift, given by the value 
ofefpr0 = TT. 

efo (1 + 6') 
V(sinhV+sin*0) (1+F^ 
Note: & is independent of the angle of 
attack. 

The ratio of the local velocity at the 
airfoil surface to the uniform stream 
velocity: 



> For airfoils whose trailing edges are not uniquely defined, it is convenient to con- 
sider the trailing edge as formed by an arcof very small radios of curvature, and to 
fair in the 1>fi carve near 6— x. 

> A convenient 20-point method is given here which in practice Is quite sufficient. 
However, any number may be used, or any interval which appears to be critical may 
be subdivided and treated separately. 



(13) 



2 



(14) c 
(15) 



0 L 



(16) F 



(17) Mp 



(18) Qu 



(19) 



■p.= £[sin(a+ aS) + sin(a+ £)]. 

The ratio of the local superstream pres- 
sure to the dynamic pressure (the term 
"superstream pressure" is used to 
designate the difference of the local 
pressure and the static pressure in the 
undisturbed uniform stream): 

«-!-($)• -d«-J,V. 

The segment of the x axis intercepted by 

the airfoil boundary. 
The lift coefficient 



0 L 



\ P cV> 



sin(a+/3) 



A point designated the "focus" of the air- 
foil. We may first define the complex 
constants c t and as 

= — f $(<t>) (cos o>+i sin o>) do> 

T 0 

^i^lf uV(o»(cos 2<t>+i sin 2<f>) d<£ 

V 0 

Then writing 



W 
we have 



1+4^+ 



and 7= -tan -1 



Ax'—Bi 2 



Then the complex coordinate of F is 

J 2 



z F = (x+iy) f 



The moment at J 1 is constant for all 
angles of attack: 
M f ~2tt p 6 s V s sin 2( 7 - 0) 

The moment coefficient referred to the 
point F: 

C3f,-^=4,r^sin2(7-j8) 
The "ideal" angle of attack: 
aj— ^ — 



where and e r denote, respectively, 
the values of e at the nose and tail; 
i.e., for 0 = 0 and 0 = *-, respectively. 
The ideal angle of attack for thin airfoils has been 
defined by Theodorsen (reference 5) as that angle for 
which the front stagnatiou point is at the leading edge. 
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At this angle of attack large velocity gradients at the 
leading edge are avoided and the profile drag is at, or 
very near, its minimum value. The definition can be 
naturally extended to actual airfoils to designate that 
angle of attack for which the front stagnation point is 
at the foremost edge of the mean-camber line. How- 
ever, as pointed out by Theodorsen, the effective mean- 
camber line of a thick airfoil actually alters with 
change of angle of attack, and the ideal angle of attack 
for a thick airfoil represents an average of a range of 
angles for which the profile drag is very near its 
minimum . 

PROCEDURE AND ACCURACY OF THE CALCULATIONS 

In order to avoid possible confusion it may be well 
to state beforehand that the term "chord" is used in 
this paper as synonymous with the segment_of the x 
axis intercepted by the airfoil. The "standard chord" 
in terms of which the airfoil is usually empirically 
denned does not, in general, coincide with this above- 
defined chord. The angle between the x axis as chosen 
and the standard chord is designated X and is listed 
in table I. (See also fig. 1.) 

In the procedure of the calculations, the axes of 
coordinates are first chosen in a definite convenient 
way, since the ease and rapidity of convergence of 
further evaluations depend considerably on the choice 
of axes (references 6 and 7). This choice may be made 
as follows: If the distance between the leading edge of 
the airfoil and the center of curvature of the leading 
edge is bisected at E (the coordinates of E are (2, 0)), 
and the same is done for the trailing edge at E' (the 
coordinates of E' are ( — 2, 0)), then the x axis should 
pass through EE' and the origin bisects the distance 
EE'. However, small variations from this particular 
choice of axis and origin do not noticeably influence 
the ease of calculation. 3 The quantities given in the 
headings of table H are directly calculated in terms of 
x and y by means of the formulas previously listed. 
The angle of attack corresponding to a given value 
of the lift coefficient may be obtained from (15), in 
which c is X^— X T , where X N and X T denote the abscis- 
sas of the leading edge.and t railing edge, respectively. 
The moment coefficient 0 VF may be obtained from 
(18), in which the constants b 2 and 7 are obtained from 
(16) by graphical integration of the sin <f>, $ cos <f>, 
■ty sin 2<j>, and f cos 2<f> curves. 

The ordinates of the airfoil are given empirically 
to hundredths of a percent of the standard chord for 
16 stations of the upper and lower surfaces respec- 
tively. The quantities x, y, \f/, and 6 are defined to 
the same degree of accuracy. The iff, 0 curve is thus a 
faired curve through 32 points and e (0) is estimated to 
be of the same order of accuracy as ^-(0). The deriva- 

• Notice, however, that we have chosen the rear stagnation point on the x axis at 
»=■*-, thus fixing the calculation by the Kutta condition. 8trlctly speaking, then, 
the z ads should be chosen In the unique manner Indicated above (this has been 
done tor the airfoils treated here), for, if another axis Is chosen, the rear stagnation 
point must be properly relocated with respect to It, causing a change In the angle of 
■rero lift with respect to the new axis. 



tives e' and $', being determined graphically, admit of 
a possible small error which, however, causes an error 
in k of probably less than 2 percent. The angle of zero 
lift, or the value of e for 8 = ir, may perhaps be in error 
as much as 15', but the influence of this possible error 
on the theoretical pressure-distribution curves for 
fixed values of 0 L is negligible. 

The numerical data for the Clark T airfoil are pre- 
sented in table H. The distribution of velocity and 
pressure for any angle of attack or at any lift coeffi- 
cient, as well as other theoretical characteristics, are 
obtained with a minimum of effort from this table. 
Similar tables for the remaining airfoils are omitted 
here for reasons of economy in printing and also be- 
cause it is not known how general the interest in them 
will be. They are available on request from the 
National Advisory Committee for Aeronautics. 

DISCUSSION 

Although the airfoils chosen in this paper are mainly 
conventional airfoils (fig. 1) and not extremely radical 
types, it is nevertheless possible to isolate some of the 
individual effects of change of shape and compare 
these with experimental results. It is believed, how- 
ever, that future experimental work on radical and 
less conventional shapes, for which the theoretical 
results are readily available (see, for example, reference 
7, p. 31), will enable the isolation and analysis of 
effects which are probably masked and unemphasized 
in conventional types. 

We may first make some general comments regard- 
ing the curves of theoretical pressure distribution given 
in the following pages. In each figure the abscissa 
represents the location of a point of the airfoil surface 
in percent chord and the ordinate gives the quantity 
y/g, the ratio of the local superstream pressure to the 
dynamic pressure g. It may be noted that negative 
values of pfa are plotted upwards. This is an arbitrary 
convention and is made because it is more readily 
associated with the upper surface of the airfoil, which 
for ordinary angles of attack is the surface of suction 
or negative pressure. In figures 2 to 21, inclusive, it 
may be noted that the points of the curves above the 
zero, or normal pressure, line represent suction; that 
is, velocities, greater than V. Positive values of ;p/g 
denote pressures greater than normal static pressure; 
i.e.. v<V. The stagnation points at which »=0 cor- 
respond to 2=1. 

Effect of compressibility. — In figure 22 there is 
shown for convenience a curve of the dynamic pressure 
g in inches of water and in pounds of force per square 
foot against velocity in miles per hour. The values 
given correspond to atmospheric conditions at 2,000 
feet altitude and 0° C- For the ordinary velocity 
range of aircraft, say from 45 to 200 miles per hour, g 
varies from about 1 to 20 inches of water. For very- 
great velocities the effect of compressibility of the air 
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1. Clark Y 




11. R.A.F.15 




Fioube L— The 20 airfoils chosen, showing axes and location of F . 



becomes significant and the potential-theory charac- 
teristics based on an assumption of incompressibility 
may be considerably altered. However, as is pointed 
out by Glauert (reference 10), the compressibility of 
air has minor influence for velocities under 0.5 the 
velocity of sound, or ordinarily about 350 miles per 
hour. However, it should be noted that at certain 
angles of attack the local velocity may be as much as 
two or more times the stream velocity. Thus, for 
very great velocities, the strong suction in the region 
of the peak pressures may introduce radical changes 
in the flow, as the compressibility properties of the 
fluid become important. This effect is associated with 
Mach's Number (u/c, where c is the velocity of sound in 
the medium), and the ordinary Reynolds Number 
alone is not a safe criterion for scale effect. The 
potential theory yet remains to be properly modified 
for the effect of compressibility. Reference 9 gives 
some experimental results of the distribution of pres- 
sure over airfoils for very high speeds. The maximum 
negative pressure (or suction) obtained in this refer- 
ence was 37 cm of mercury. 

Pressure gradients. — From the concept of the ideal 
angle of attack we are led to expect that the thinner the 
airfoil at the leading edge the greater the velocities 
near the leading edge for angles different from a r . This 
expectation is confirmed by the large negative values 
of fig attained by the RA.F. 15 and the N.A.O.A. 
0010, 0012, and 2409 airfoils. In particular, the pres- 
sure on the N1.C.A. 0010 reaches — llg for lift 
coefficient Cx=1.5, whereas the somewhat thicker 
NA.C1. 0012 reaches — 7q at the same lift coefficient. 
In practice, the value of C Laax for the N.A.C.A. 0012 
is somewhat greater than that for the N.A.C.A. 0010. 
Results of force tests of the airfoils treated in this 
paper are presented in references 11, 12, and 13. 

The large gradient of pressure behind the negative 
peak pressure is very significant for the breakdown of 
potential flow. The deceleration of fluid becomes so 
rapid that fluid is piled up at the trailing edge and the 
flow no longer separates precisely at the trailing edge 
but breaks off along the upper surface. The flow along 
the lower surface undergoes but little change at high 
angles of attack except that, after the breakdown of 
potential flow has occurred, the pressure at the trailing 
edge may be somewhat negative instead of positive. 

The change from the front stagnation point to maxi- 
mum velocity occurs within a very small space interval, 
and all indications are that frictional losses are practi- 
cally negligible while the fluid is accelerating, as com- 
pared with losses when deceleration occurs. The 
fluid follows the surface boundary more easily in the 
change from pressure energy to kinetic energy than 
vice versa. This fact is also abundantly confirmed by 
experiments with nozzles and venturi tubes. For 
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Figure 2.— Theoretical pressure distribution tor the Clark Y airfoIL 




Figure 3.— Theoretical pressure distribution for the N.A.O.A.-M6 airfoIL 




Figure 4.— Theoretical pressure distribution tor the TJ.S.A. 27 airfoIL 




Figure 6.— Theoretical pressure distribution for the TT.8.A. 36B airfoIL 
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Fiocre 6.— Theoretical pressure distribution for the N 22 airfoil. 
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Fiocee 7.— Theoretical pressnre distribution for the GBttlngen 387 airfoil. 
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Figure a— Theoretical pressnre distribution for the GBttlngen 398 airfoil. 
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Figure 9.— Theoretical pressnre distribution for the N.A.C.A.-CYH airfoil. 
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Fiouee 11.— Theoretical pressure distribution for the Boeing lOS-A airfoil. 




Figtjee 12.— Theoretical pressure distribution for the R.A.F. 15 airfoil. 
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Fiotoe 13.— Theoretical pressure distribution for the R.A.P. 19 alrfoD. 
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FiauttB 14.— Theoretical pressure distribution for the N.A.O.A. 0010 airfoil. 
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Figure 17.— Theoretical pressure distribution for the N.A.C.A. 3409 airfoil. 




Fionas 18.— Theoretical pressure distribution for the N.A.O.A. 2412 airfoil. 
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airfoils used in practice without auxiliary devices, a 
pressure ratio fig of about —3 or —4 is the m fori mum 
attained. It may be observed here that since the 
leading edge of the airfoil and the upper surface near 
the leading edge experience very large pressure 
gradients, they are critical regions, to be especially kept 
free from unnecessary protuberances and roughness. 

It may be noted in the figures that very large 
pressure gradients exist in the region near the rear 



Some effects of camber.— It may be observed that a 
property common to all airfoils is that the negative 
values of fig mount rapidly near the leading edge on the 
upper surface after Ox, =1.00 is exceeded. However, 
for the highly cambered airfoils, as the R.A.F. 19 and 
N.A.C.A. 6512, it may be noted that even for 0 L = 1.50 
the f/g negative peak is but slightly above — 2. These 
airfoils are high-lift airfoils, and the lift is well distri- 
buted over the whole chord. However, they have 
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Figure 19.— Theoretical pressure distribution for trie N.A.O.A. 2416 airfoil. 
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Figure 20.— Theoretical pre s ume distribution for the N.A.O.A. 4412 airfoil. 
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Figure 2L— Theoretical pressure distribution for the N.A.C.A. 6612 airfoil. 



stagnation point, within the shadow of the 100 percent 
chord line. The rapid deceleration of fluid as thus 
shown to exist theoretically at the trailing edge most 
probably does not occur to any such extent in practice. 
The flow probably recombines at the trailing edge at 
velocities not much below normal, and the pressure 
curves are rounded off at a small positive pressure as 
shown in the figures. There is in this fact no essential 
violation of the Kutta condition for fixing the circula- 
tion, the primary purpose of which is only to avoid 
infinite velocities at the trailing edge. 



usually unfavorable pitching moments and rather wide 
travel of the center of pressure, as evidenced by the 
theoretical moment coefficients at zero lift, which are 
respectively —0.210 and —0.185. 

Further effects of camber may be observed in figures 
15, 18, 20, and 21 for the N.A.O.A. 0012, 2412, 4412, 
and 6512 airfoils, where in every case the maximum 
thickness is 12 percent of the chord and the maximum 
mean cambers are, respectively, 0, 2, 4, and 6 percent 
nf the chord. The theoretical moment coefficients 
CWare, respectively, 0, -0.055, -0.110, and -0.185. 
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The ideal angle, also, increases with camber and, hence, 
the optimum lift coefficient increases in general with 
camber. We may observe that the high-cambered air- 
foils, and more especially thin high-cambered airfoils, 
are not efficient at low values of the lif t coefficient, as is 
evidenced by the high negative peaks in the pressure 
distribution for zero lift. In fact, the flow around 
many high-cambered airfoils (for example, the R1.F. 
19) is known to burble on the under surface at low 
lift coefficients. Indeed, large gradients of decelera- 
tion of fluid may everywhere be associated with de- 
creased efficiency. The bringing of pressure gradients 
and profile resistance into a precise correlation is a 
significant problem for further investigation. A 
uniform gradual change from velocity to pressure, as in 
figure 2 for the Clark T airfoil at C L =0.5, gives prob- 
ably the optimum lift distribution and occurs very 
nearly at the ideal angle of attack. 

The experimental curves of lift coefficient against 
angle of attack for high-cambered airfoils like the 
R.A.F. 19 and N.A.C.A. 6512 are well rounded near 
maximum lift, whereas for airfoils like the N.A.C.A. 
0010 and R.A.F. 15 they are likely to be sharp and 
jagged (reference 12). The former airfoils lose their 
lift rather gradually after maximum lift is attained, 
while for the latter airfoils this effect is likely to occur 
suddenly. 

Effects of thickness. — In figures 17, 18, and 19 we 
may note some effects of the airfoil thickness. The 
N.A.C.A. 2409, 2412, and 2415 airfoils have a common 
mean-camber line and maximum thicknesses, respec- 
tively, of 9, 12, and 15 percent of the chord. For the 
N.A.C.A. 2415 it may be noted that the pressure on 
the under surface is generally less positive than for 
the 2412 and 2409. Also, we may observe that at lift 
coefficients C L = 1.00 and Ok = 1.50 the down gradient 
of pressure along the first 15 percent of the chord is 
greatest for the 2409, while for the rest of the chord it 
is greatest for the 2415, indicating that an optimum 
effect for thickness lies perhaps between the extremes 
listed. The theoretical slope of the lift curve in- 
creases somewhat with thickness, and for the above 
airfoils has values of 6.75, 6.90, and 7.10, respectively. 
An experimental result that merits closer investigation 
is the fact that after a maximum thickness of about 12 
dO 

percent is attained, decreases somewhat with fur- 
ther increase in thickness (reference 12). A partial 
explanation lies in the fact that, in general, a thicker 
boundary layer exists on thick airfoils, decreasing their 
aerodynamic efficiency. It has, indeed, puzzled some 
observers employing various schemes for removing the 
boundary layer that the experimental slope of the lift 
curve for infinite aspect ratio often exceeds 2a-, which 
is the value given by the approximate thin-airfoil 
theory. The average theoretical slope of the curves 
of lift against angle of attack for the airfoils listed in 



this report is approximately 6.90 oi about 1.10X2tt. 
This value is about 15 to 20 percent greater than the 
experimental value of the slope of the lift curve for 
infinite aspect ratio, and indicates that the airfoil is in 
general from about 80 to 85 percent efficient with 
regard to lift. 

Moment properties. — The theoretical moment 
M F is, in most of the cases studied in this paper, from 
about 10 to 20 percent greater than the experimental 
value of the moment for zero lift. (See references 12 
and 13.) On the basis of (16), the position of the focus 
F, at which the theoretical moment is constant for all 
angles of attack, was calculated and is shown in the 
figures of the airfoils. (See fig. 1 and also table I.) 
In every case its abscissa is very nearly at 25 percent 
of the chord from the leading edge, the maximum for 




V, velocity in miles per hour 
Fioubz IX— Velocity V against dynamlo pressure ff. 



the airfoils considered being near 27 percent and the 
minimum near 24 percent of the chord line. How- 
ever, it is important to note that, in general, the ordi- 
nate of F does not fall on the chord fine but is usually 
located at a small distance from it. 

The tendency for more constant center-of-pressure 
properties may be observed in figure 4 for the U.S.A. 
27, where the lower-surface pressure curve has a small 
inflection or bend. For the N.A.C.A.-M6 (fig. 3), 
which theoretically has practically zero travel of the 
center of pressure, the double bend on the pressure 
curves at zero lift may be observed. Alternate re- 
gions of suction and pressure thus exist on each surface. 
The N.A.C.A.-CYH shows the same tendencies to a 
lesser degree. The double bend in the pressure- 
distribution curves is probably common to most 
reflexed airfoils. 
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Experimental results on the angles of zero lift show 
considerable discrepancies and indicate a change with 
the Reynolds Number. In general, the consistent 
experimental result is obtained that the angle of zero 
lift increases (algebraically) with increase in the Rey- 
nolds Number. For this reason we may only indicate 
a comparison with experimental values. The values 
listed in table I of this paper consistently fall within 
the range of values obtained by experiment and seem 
to agree more closely with experimental results obtained 
at moderate Reynolds Numbers (about 2 X10 6 ) than 
with those taken at very large Reynolds Numbers. 
This fact should be investigated further. 

CONCLUSIONS 

The preceding discussion shows that, to a large 
extent, the general properties and characteristics of 
airfoils, such as effects due to camber, thickness, or 
change of shape, may be accounted for by the theoreti- 
cal pressure-distribution curves. The theoretical pres- 
sure-distribution curves at definite lift coefficients may 
be safely used for structural-load considerations. A 
correlation between pressure gradients and profile 
resistance, affecting the efficiency of the airfoil, has 
been qualitatively indicated. Further theoretical and 
experimental investigations should be concerned with 
the significance of the differences between theory and 
experiment. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 2, 19SS. 
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